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Abstract 

Materials with strong spin-orbit coupling have attracted attention following the 
prediction and subsequent discovery of strong two- and three-dimensional topological 
insulators in which a topological property of the bulk band structure of an insu- 
lator results in metallic surface states with Dirac-like dispersion. Here we report 
the discovery of Dirac-like surface states in the perovskite superconductor [T^jTlTea 
(TlsTes) and its non-superconducting tin-doped derivative, [Tl4](Tlo.4Sno.6)Te3, as ob- 
served by angle-resolved photoemission spectroscopy (ARPES). Density functional 
theory (DFT) calculations predict a single spin-orbit driven band parity inversion at 
the Z point above the Fermi level of Tl5Te3, suggesting the surface states are protected 
by Z2 topology. Calculations on [Tl4]SnTe3 show no parity inversions, implying that 
a topological transition from non-trivial to trivial must occur upon doping with tin, 
i.e., [Tl4](Tli^Snj:)Te3. Thus [Tli]MTe3 perovskites are a possible new, non-trigonal 
class of Z2 topological compounds. Additionally, as TlsTea is a stoichiometric bulk 
superconductor, these perovskites are ideal materials in which to study the interplay 
between surface states and bulk superconductivity. 
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Charge carriers in Dirac-like bands have a hnear energy-momentum relationship described 
by the Dirac equation; they offer the unique opportunity to investigate the intersection of 
special relativity and quantum mechanics. However, materials exhibiting Dirac surface states 
are relatively rare, the most notable examples being graphene flj, Z2 topological insulators 
[2], and topological crystalline insulators (TCIs) [3J. Experimentally known topological 
insulators are based on a limited range of materials and structure types: strained HgTe, 
Bi^Sbi_^, (Bi/Sb)2(Se/Te)3, Tl[Bi/Sb](S/Se/Te)2, and 3114^3X9 gHH] with many others 
theoretically predicted P4T^. Materials which combine topological surface states with su- 
perconductivity have the ability to test the predictions of new states of matter and exotic 
quasiparticles, such as Majorana fermions, and open the field for new technologies such as 
quantum computing [T3H16]. Producing such materials has been a challenge: the known 
methods of combining these two properties are proximity effects and doping (Cu-doping of 
Bi2Se3 and In-doping of Sui _^In^Te) p!71 - [T9] . which limits studies to small samples and/or 
materials with low volume fraction superconductivity. 

Here we report the discovery of Dirac-like surface states in the tetragonal perovskites 
[Tl4](Tli_2.Sni.)Te3, as observed by angle-resolved photoemission spectroscopy (ARPES). 
The surface states are predicted by density functional theory (DFT) to be Z2 topologically 
protected due to a spin-orbit-driven band parity inversion at the Z point. This is the 
first observation of Dirac-like surface states in a material with the perovskite structure 
and represents a possible new, non-trigonal class of Z2 topological compounds. Moreover, 
we show that ~100% volume fraction superconductivity is possible in Tl5Te3, making it 
an ideal compound in which to explore the coupling between bulk superconductivity and 
Dirac-like surface states, i.e., a "topological superconductor." 

RESULTS 

Angle-Resolved Photoemission Spectroscopy 

The [TUjMTeg (M = Sn, Pb, Bi, Sb, La, Nd, and Mo [20H23]) family of compounds 
has perovskite structure (AiJOs) described by Glazer notation a°a°c~ with body-centered 
tetragonal symmetry pH [25]. The cavities of a three-dimensional network of corner-sharing 
MTee octahedra are occupied by interconnected TI4 tetrahedra [Fig. [l]^a)]. 
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In Fig. [2]^a) we show ARPES data on [Tl4](Tlo.4Sno.6)Te3. The most prominent feature is 
an intense, hnearly dispersing, Dirac-hke band centered at a momentum of (k^, = 0,ky = 0) 
of the two-dimensional surface Brillouin zone (2DBZ). We observe this feature at the center 
of numerous Brillouin zones, and only a single such metallic feature in each zone. We assign 
this Dirac-like band to be a surface state based on a number of different reasons. First, 
the values of the extracted Fermi velocity [4.6(2) eV A] and the Fermi wavevector [0.04(1) 
A~^] are insensitive to incident photon energy indicating that they are non-dispersive in 
kz [Fig. [2]^c)], even though the bulk electronic structure is highly three-dimensional. This 
is a hallmark of a surface state. Furthermore, this feature is present and always centered 
at (0,0) of the 2DBZ, and not the center of the three-dimensional bulk Brillouin zone, as 
would be expected from a bulk-derived state. Additionally, we also observe features near the 
Fermi energy {Ep) that do vary considerably with the incident photon energy, and which 
we ascribe to bulk states, consistent with DFT calculations (see below). 

In Fig. |3| we show ARPES data from TlsTes. Much like in [Tl4](Tlo,4Sno.6)Te3, a single, 
Dirac-like surface state is once again visible at (0,0) while no other Dirac-like features are 
observed anywhere else in the 2DBZ. As TlsTes is metallic, containing fewer electrons than 
[Tl4](Tlo.4Sno.6)Te3, only the bottom section of the Dirac cone is visible, and the extrapolated 
Dirac point sits 200 meV above E^?. The Fermi velocities [3.0(2) eV A] and Fermi wavevectors 
[0.07(1) A~^] observed for the Dirac-like band in TlsTcs are largely independent of photon 
energy {hu = 21.2 and 40.8 eV). 

The similarities between surface states observed both in [Tl4](Tlo.4Sno.6)Te3 and Tl5Te3 
and the single Dirac cone observed per Brillouin zone suggest that [Tl4](Tli_2:Sn2:)Te3 host 
topologically protected surface states. Further indirect evidence for their topological nature 
is that we find no Rashba-type spin splitting of the surface state despite the strong spin-orbit 
coupling of these materials, suggesting that we have just a single spin non-degenerate Fermi 
surface between neighboring time-reversal invariant points, the hallmark of a TI. 

Density Functional Theory Calculations 

However, there are several ways in which the [Tl4]MTe3 structure can result in topological 
Dirac-like surface states. First, the mirror planes which topologically protect the Dirac states 
in rock salt TCIs are also present in the basic perovskite structure. The appropriate band 
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inversions at points protected by mirror symmetry similarly may yield a TCI state in body- 
centered tetragonal symmetry. Second, the body-centered tetragonal crystal class has four 
sets of time-reversal invariant points: one F, four A^, two X, and one Z. If spin-orbit 
coupling produces band parity inversions at the F and Z points, but not or X, surface 
states protected as in a weak topological insulator result. Alternately, a band inversion due 
to spin-orbit coupling at the T or Z point, but not both, produces strong Z2 topologically 
protected surface states. 

Calculations for TlsTcs, shown in Fig. [4], show spin-orbit coupling induces a single band 
inversion approximately 0.5 eV above the Fermi level between valence and conduction band 
states at the Z point. This inversion is between opposite-parity states, the conduction band 
involving states from the [TI4] units on the A site and the valence band involving mainly 
states from apical Te atoms and B site Tl ions in the corner- sharing perovskite network. 
If electron-doped to the gap at Z opened by spin-orbit coupling, this system would be a 
topological semi-metal due to high- lying valence band states at the F point. The resulting 
situation is thus similar to that found in pure antimony [S]. Furthermore, we have calculated 
the Z2 invariant to be -1 by taking product of the parity of all the bands at the F and Z 
points that connect to the occupied valence bands, including two bands above the Fermi 
level at the F point. Therefore, DFT predicts TlsTes to have surface states protected by a 
non-trivial Z2 topology. 

In contrast, our calculations shown in Fig. [s] demonstrate that [Tl4]SnTe3, the opposite 
end member of the [Tl4](Tli_a.Sna.)Te3 solid solution, has trivial Z2 topology: although a 
possible band inversion occurs near the F point, the involved states are of the same parity. 
Further, an approximately 0.7 eV band gap opens up at the Z point where the band crossing 
is observed in TlsTes, preventing spin-orbit coupling from inverting band parity at that point. 
These calculations imply that a transition from Z2 non-trivial (at x =0, or TlsTca) to trivial 
(at X =1, or [Tl4]SnTe3) topology must exist at some doping level x in [Tl4](Tli_^.Sna,)Te3 
compounds. The situation in the [Tl4](Tli_a;Sna;)Te3 peroskites is thus similar to what is 
found for the Bii_i.Sbx. solid solution [5]. 

To understand the driving force behind the topological transition, we investigated the 
origins of the large gap opening at the Z point when Sn is substituted for Tl. There are 
two major structural changes when this doping is performed: a sizable change in lattice pa- 
rameters, particularly the increase in the c axis relative to a, and cation mixing. In order to 
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isolate the effect of lattice parameters, the band structure of TlsTcs was calculated with the 
same lattice constants as [T^jSnTes imposed, and it was found that TlsTcs would then be 
topologically trivial, albeit with a smaller bandgap at the Z point. Further emphasizing the 
importance of the lattice parameters on the topological transition is the fact that the width 
of the valence band involved in the band inversion that has large apical Te weight depends 
sensitively on the lattice constant c of the material due to slight B site Tl contribution. This 
therefore demonstrates the importance of lattice parameters in determining the Z2 topolog- 
ical invariant in the [Tl4](Tli_j.Snj.)Te3 compounds. Preliminary calculations (not shown) 
suggest that the effect of cation ordering is also significant, but detailed studies beyond the 
scope of this work are needed to elucidate those details and to determine the critical value 
of X required for the transition to a different topological phase in [Tl4](Tli_2.Sn2.)Te3 com- 
pounds. The sensitivity to structural parameters, particularly the c/a ratio, suggest that 
temperature- and/or strain-driven topological transitions may be possible in this material 
family. Regardless of these details, DFT results are consistent with the ARPES detection 
of a Dirac cone in both the TlsTca and [Tl4](Tlo.4Sno.6)Te3 compounds and indicate that Z2 
topology protects surface states in this materials family. 

Superconducting Properties 

The observation of Dirac-like surface states in TlsTes, likely Z2 topologically protected, 
makes it a candidate topological superconductor. Superconductivity was not observed in 
[Tl4](Tlo.4Sno.6)Te3 above T = 1.8 K. Because Tl5Te3 was previously proposed to be a two- 
band superconductor [261 EZ] , we performed a comprehensive specific heat and magnetization 
study of the superconducting properties. Our specific heat data collected on a polycrystalline 
chunk under zero magnetic field reveal a clear A-anomaly. An equal-entropy construction, 
shown in Fig. |6](a), indicates bulk superconductivity with a Tc = 2.40(1) K. The ratio 
ACez/(7Tc) = 1.63 is close to, but slightly greater, than the weak coupling Bardeen-Cooper- 
Schrieffer (BCS) value of 1.43 [28]. The normal state heat capacity, measured under H = 5 
kOe, is well described by C„ = 7T + PsT^ + P^T^, where the term linear in T represents the 
electronic contribution and the higher order terms represent the lattice contribution [2S]. 
The fit yields an electronic Sommerfeld coefficient of 7 = 10.6(1) mJ mol f.u.~^ K~^. The 
coefficients P3 and /Ss of the lattice contribution are 3.90(5) mJ mol f.u.~^ and 0.740(6) 
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mJ mol f.u. ^ K ^, respectively, and the /Ss corresponds to a Debye temperature of 6d = 
159(2) K [23. 

Low-temperature specific heat measurements collected under H = Oe and 5 kOe show 
the latter field is sufficient to quench superconductivity. The data collected under a. H = 
5 kOe field below T = 1.5 K, temperatures at which any complex phonon structure should 
be negligible, were fit to Cei = 7T + /SsT^, where 7 = 9.34(4) mJ mol f.u.~^ and /Ss 
= 6.82(5) mJ mol f.u.~^ K""^ [29]. Subtracting this low-temperature lattice specific heat fit 
gives the electronic specific heat, which is well described by the BCS model of the electronic 
specific heat of a single-gap s-wave superconductor Cei/{jTc) = Aexp[(-BTc)/T] with A = 
11.4(1.1) and B = 1.61(2) [Fig. gb)] [2HI|29]. According to BCS theory, B = A ^bTc, and 
thus a value of B = 1.61(2) predicts a superconducting gap of A = 0.333(6) meV. This 
value, based on the electronic specific heat data, agrees with the superconducting gap of 
A = 0.365(2) meV predicted using the critical temperature of Tc = 2.40(1) K in the BCS 
relation 2A/ (k^Tc) = 3.53. 

Isothermal DC magnetization measurements were collected on an approximately cubic 
100 mg single crystal of TlsTcs with (001) and (110) faces, allowing the orientation depen- 
dence to be examined. At T = 1.8 K, a full superconducting volume fraction is observed 
from both crystal directions under small applied magnetic field, and the upper critical fields, 
Hc2,a6(l-8 K) = 480(30) Oe and Hc2,c(1.8 K) = 500(50) Oe are within error. In contrast, the 
lower critical fields are found to differ by a factor of greater than three, as shown in Fig. 
gc): H*hi'^.8 K) = 18(1) Oe and H*{1.8 K) = 5(1) Oe, where H*{1.8 K) is defined as the 
magnetic field, parallel to the i direction, at which the magnetic response deviates from a 
linear field-dependence. We call this H* rather than the strictly defined lower critical field 
Hci^i because the presence of surface states may produce spurious effects in the measured 
magnetic response. Note this observed threefold anisotropy is merely a lower bound of the 
extent of the anisotropy present in the bulk superconducting state: any intrinsic material 
defect of the crystal or minor experimental misalignment would result in contribution from 
the perpendicular direction. In short, TlsTca appears to be a single-gap BCS supercon- 
ductor. It is thus an ideal material in which to probe whether topological surface states 
participate in superconductivity, or how these gap less states affect, and are affected by, the 
bulk superconducting gap. 
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DISCUSSION 



Charge carriers in Dirac-like bands have a hnear energy-momentum relationship described 
by the Dirac equation; they offer the unique opportunity to investigate the intersection of 
special relativity and quantum mechanics. The Dirac cones of graphene are even-numbered 
and topologically trivial. In contrast, the surface states of Z2 topological insulators and 
TCIs are protected topologically by time-reversal and mirror symmetry, respectively. Strong 
Z2 topological insulators must have an odd number of band inversions at time-reversal 
invariant points to ensure protected surface states. Electron transport on the surface of a 
Z2 topological insulator is novel because the spin and momentum of surface electrons are 
"locked" together. For this reason, topological insulators have a wide variety of potential 
applications. 

Our results suggest a number of characteristics which make the [Tl4]MTe3 family an in- 
triguing alternative to most other topological insulator families. First, as shown in Fig. |2](d). 




the truncation of the Dirac cone indicates that [Tl4](Tlo.4Sno.6)Te3 is p-tjpe doped, in con- 
trast to materials such as Bi2Se3, which are typically ra-type in the bulk. Second, instead of 
the six-fold warping of the Dirac cone away from the Dirac point which is consistent with 
the trigonal structure of most topological insulators, the Dirac cone in [Tl4](Tlo.4Sno.6)Te3 
exhibits a four-fold warping consistent with its tetragonal symmetry. At higher binding 
energies, the surface state appears to intersect and begin to hybridize with the bulk bands 
[such as around 0.6 eV in Fig. [2|^a)], forming a surface resonance with some k^ dependence 
in the Fermi velocity. However, at energies near Ep within the bulk bandgap, this band 
exhibits a slightly higher Fermi velocity shown in Fig. ^c). Third, DFT analysis shows that 
the topological transition from trivial to non-trivial must occur in [Tl4](Tli_a;Sna;)Te3 as x 
decreases. The sensitivity of this transition to structural parameters, particularly the c/a 
ratio, suggesting that temperature-, pressure-, or stain-induced topological transition may 
be possible in this material family. 

The coexistence of topologically protected Dirac-like surface states and superconductivity 
in TlsTes make it a candidate "topological superconductor." Topological superconductivity 
is highly sought after due to the prediction of massless Majorana fermions, which have 
great potential to advance technology, specifically in the area of quantum computing, and 
to improve our knowledge of fundamental physics. Although superconductivity has previ- 
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ously been achieved in topological systems via proximity effects and by doping, experiments 
investigating the phenomenon are hindered by several factors. For example, while super- 
conductivity has previously been realized in Cuj^Bi2Se3 by intercalating copper into the Van 
der Waals gaps between layers, a small superconducting volume fraction is observed in this 
material [18]. Further, the substitution of copper for bismuth (Bi2_xCuj;Se3) competes with 
copper intercalation (Cua;Bi2Se3), thereby limiting the attainable Cu dopant concentration 
and producing inhomogeneous samples [18] . Superconductivity in topological insulators has 
also been successfully induced in Bi2(Se/Te)3 by the proximity effect, but studies are then 
limited to thin sample regions[I7j. Sni_j.In^.Te and LaPtBi are also candidate materials, 
but the former is a TCI and the latter lacks a center of inversion, complicating efforts to 
understand them [THl US] . 

In contrast to these materials, the [T^jMTes family provides an ideal system to investi- 
gate the effects of topologically protected surface states on superconductivity because TlsTea 
is an intrinsic bulk superconductor that is straightforward to synthesize and computation- 
ally model. Tl5Te3's bulk superconductivity eliminates the experimental obstacles faced 
in doped or proximity-induced topological superconductors, such as low superconducting 
volume fraction, sample inhomogeneity, and size limitations. Further, the wide range of 
cations reported to substitute for the M site, as well as the possibility of substituting other 
chalcogens or halogens for tellerium, presents a versatile range of experimentally accessible 
doping levels and resultant properties to examine. 

The [Tl4]MTe3 family has an entirely different crystal structure than the known Z2 topo- 
logical insulators, which all have trigonal symmtry in the active units. Because symmetry 
determines the type and number of time-reversal invariant momentum points in the Bril- 
liouin zone of a material as well as its mirror planes, the structure type of a material is an 
important factor in whether Dirac surface states exist and whether they are topologically 
non-trivial. Our finding of topological surface states in perovskites therefore has impor- 
tant materials implications. Other members of the [Tl4]AfTe3 family likely exhibit similar 
topological properties in conjunction with other emergent phenomena including supercon- 
ductivity. More generally, our results open the door to finding topologically protected states 
in other members of the largest known material family - perovskites. 
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METHODS 



Sample Preparation 

Polycrystalline TlsTes samples were prepared by heating elemental Tl (Strem Chemicals, 
99.9 %, warning: toxic, carcinogenic) and Te (Alfa Aesar, 99.999+ %) in a vacuum-sealed 
silica ampoule to 550°C, holding for 24 h, followed by slow cooling (5°C/h). Elemental 
thallium was stored and handled in an inert atmosphere to prevent oxidization. Single 
crystals were prepared using a modified Bridgman method in an optical floating-zone furnace 
(Crystal Systems Inc.) using 2.5% excess Te as a flux. Tl4(Tli_j,.Sn^.)Te3 samples were 
prepared by heating elemental Tl, Te, and Sn (NOAH Technologies, 99.9%) in a vacuum- 
sealed silica ampuole to 540°C, holding for 24 h, followed by slow cooling (2°C/h). Single 
crystals of Tl4(Tlo.4Sno.6)Te3 suitable for ARPES measurements were obtained from the 
ingot after slow cooling. 

Physical Properties Measurements 

Physical property measurements were performed on a Quantum Design, Inc. Physical 
Property Measurement System. Low-temperature heat capacity data from T = 0.15 K to 4 
K were collected on a polycrystalline chunk with the aid of a dilution refrigerator. ARPES 
measurements were performed at the Synchrotron Radiation Center in Wisconsin at the 
U3 beamline with a VG Scienta R4000 electron analyzer at a base temperature of 30 K 
and a pressure of 8 x 10^^^ torr as well as a lab-based ARPES system with a VG Scienta 
R4000 analyzer and VUV5000 helium plasma discharge lamp {hu = 21.2 and 40.8 eV) and 
monochromator at a temperature of 10 K and a base pressure of 6 x 10^^^ torr. Samples 
were cleaved in situ along the c axis in ultra-high vacuum in order to expose atomically 
clean surfaces. Due to the three-dimensionally bonded structure, cleaved samples did not 
yield entirely fiat surfaces, but still possessed regions with reasonably large fiat facets. 

Density Function Theory Calculations 

First principles (DFT) calculations were performed with full-potential linear augmented 
wave formalism as implemented in WIEN2k [30]. Exchange-correlation energies are calcu- 
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lated with Perdew-Burke-Erzenhof (PBE) functional [31]. An 8x8x8 unshifted k-point grid 
is used in the Brillouin Zone of the primitive body-centered tetragonal cell. Results re- 
ported were double-checked by repeating calculations in Projector Augmented Wave (PAW) 
formalism, as implemented in VASP, and no disagreements were found [32] . 
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FIG. 1: (Color online) Structural model of the perovskite framework of [Tl4]MTe3 (M = 

metal). TI4 tetrahedra occupy the A site of a conventional perovskite ABO3, and the alternating 
rotation of AfTeg octahedra about the c axis is described as a^a^c~ in Glazer notation with a 
rotation angle of 20.1° for TlgTes [Ml ESI [33] . 
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FIG. 2: (Color online) Dirac-like surface state in [Tl4](Tlo.4Sno.6)Te3. (a) ARPES spectrum 
taken along the (0,0) - (vTjTt) shows the electronic states over a wide energy range, including 
bulk states, and a linearly dispersing surface state which evolves into a surface resonance upon 
intersecting the bulk states, (b) Expanded view of the Dirac-like surface state within the projected 
bulk band gap. (c) Fits to the dispersion at various photon energies, revealing no substantial 
dispersion along k^. (d) Momentum-resolved intensity maps plotted at different binding energies, 
which suggest the bottom portion of a single Dirac cone centered at the zone center, consistent 
with a topological surface state, with a square warping at higher binding energies, reflecting the 
tetragonal crystal symmetry. 
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FIG. 3: (Color online) Dirac-like surface states in the topological superconductor TlsTea. 

(a) ARPES spectrum taken with He la {hv = 21.2 eV) photons along the (0,0) - (7r,7r) direction 
showing the coexistence of prominent bulk states and a linearly dispersing surface state (yellow 
dashed lines) near Fip. (b) Second derivative plot of the measured experimental data in a, where 
the surface state is more clearly resolved, (c) MDC fits for the data taken with He la and He 
Ha (hi/ = 21.2 eV and 40.8 eV) photons, showing no substantial dependence on photon energy, 
(d) Momentum-resolved intensity maps at a series of binding energies which show the single Dirac 
cone feature centered at the zone center. 
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FIG. 4: (Color online) Spin-orbit driven band inversion in TlsTes.The band structure of 
TlsTea, calculated from first principles DFT, (a) without and with (b) spin-orbit interactions. 
Shading denotes bands arising from valence (light blue) and conduction (light red) band states, 
(c) The spin-orbit interaction produces a single band parity inversion and opens a small gap at the 
Z point, indicating that TlsTea would be a topological semi-metal if doped to the gap. (d) The 
Brillouin zone of the body-centered tetragonal unit cell for comparison with the ARPES results. 
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FIG. 5: (Color online) Topologically trivial nature of [Tl4]SnTe3. The band structure of 
TUSnTes, calculated from first principles DFT, is shown (a) without and with (b) spin-orbit 
interactions. Shading denotes bands arising from valence (light blue) and conduction (light red) 
band states. Unlike the TlsTes case, the spin-orbit interaction does not produce a band parity 
inversion at the Z point, implying that [TUJSnTes is topologically trivial semi- metal. The band 
inversion near the F point is between bands of the same parity and does not affect topological class. 
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FIG. 6: (Color online) Superconducting properties of topological TlsTea. (a) Specific heat 
data collected under zero field (purple triangles) show a clear A-anomaly at = 2.40(1) K, while 
those collected under an applied 5 kOe field (gray squares) fit well to C„ = 7T + P^T^ + P^T^- 
Inset: The electronic specific heat as a function of temperature is shown along with an equal en- 
tropy construction illustrating bulk superconductivity with Tc = 2.40(1) K. (b) Low-temperature 
electronic specific heat behaves in accordance with single-gap BCS theory, (c) Isothermal magne- 
tization data collected on an oriented single crystal at T = 1.8 K with the field applied parallel to 
the [001] direction (red diamonds) and [110] direction (blue circles) reveal a striking anisotropy of 
H*. The slope of the data in the linear region below H* indicate a 96% superconducting volume 
fraction. A demagnetization factor N = 1/3 has been applied to account for the approximately 
cubic shape of the 100 mg single crystal sample, (d) The field-temperature phase diagram shows 
Hc2 as determined by specific heat (purple triangles) and orientation-dependent isothermal DC 
magnetization (blue circles and red diamonds), along with a two-fluid model fit to the data. The 
orientation anisotropy present in !{*■ is not observed for Hc2- 
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